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Abstract 

We investigate the influence of the full vacuum polarization and vertex func- 
tion on the fermion propagator, using the coupled Dyson-Schwinger equations 
for the photon and fermion propagator. We show that, within a range of ver- 
tex functions, the general behavior of the fermion propagator does not depend 
on the exact details of the vertex, both in the massless and in the massive 
phase. Independent of the precise vertex function, there is a critical number of 
fermion flavors for dynamical mass generation in (2 + l)-dimensional QED. A 
consistent treatment of the vacuum polarization is essential for these results. 
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I. INTRODUCTION 



Quantum electrodynamics in 2 space- and 1 time-dimension, QED3, has several interest- 
ing features. It exhibits dynamical mass generation [|T| p~0|| and confinement pd]-p^|, similar 
to QCD. Furthermore, it is superrenormalizable, so it does not suffer from the ultraviolet 
divergences which are present in QED4. The coupling constant is dimensionful, and pro- 
vides us with a mass scale, even if we consider massless fermions. This energy scale plays 
the role of the QCD scale Aqcd, in the sense that it sets the scale for confinement and dy- 
namical mass generation. Thus it is a very interesting model to study these nonperturbative 
phenomena. 

QED3 also has some applications in condensed matter physics, where it can be regarded 
as an effective theory for more realistic microscopic models [JT^JTSI. Especially, since the 



discovery of high-T c superconductivity and the fractional quantum Hall effect, these kinds 
of models have attracted more attention. 

In this paper, we consider QED3 with N fermion flavors of four- component spinors. Such 
a model of QED3 is chirally symmetric in the absence of a bare fermion mass term, m O '0'0, 
in contrast to the (2+l)-dimensional gauge theory with two-component fermions, where we 
cannot define chiral symmetry |l]||,|n|. Similar to the four-dimensional case [[17|], the chiral 



symmetry of QED3 may be broken spontaneously due to the dynamical generation of a 
fermion mass. The question of whether or not chiral symmetry is broken for all values of N, 
the number of fermion flavors, is very interesting. 

We address this question by analyzing the behavior of the full fermion propagator nonper- 
turbatively, using its Dyson-Schwinger equation. However, this equation cannot be solved 
without truncating the infinite set of Dyson-Schwinger equations, since it involves the full 
photon propagator and the full vertex function. Different ways of truncating this equation 
give rise to different results concerning the question of whether there is a critical number of 
fermion flavors for dynamical mass generation 0|||6|-[|. In this paper we try to resolve the 
controversy. 

Using the bare vertex approximation and the one-loop vacuum polarization, Appelquist 
et al. have shown that there is a finite critical number of flavors N c , above which the chiral 
symmetry is restored. They found in Landau gauge a critical value of N c = 32 /n 2 ~ 3.24. 
This approach is based on a l/N expansion, and including the next-to- leading-order terms 
||, it was found that the critical number changes to N c = 128/37T 2 ~ 4.32. 

However, such a simple treatment was criticized by Pennington et al. ||, since the 
effects of the wavefunction renormalization were not taken into account. The problem is 
that formally the wavefunction renormalization is of order 1 + 0(1/N), but in the infrared 
region, this wavefunction renormalization tends to vanish. Because of this non-uniformity 
in the l/N expansion, an approximation based on an expansion in l/N might not be very 
reliable. Taking into account wavefunction renormalization, Pennington et al. found chiral 
symmetry breaking for all numbers of fermion flavors. However, as we show, their approach 
also has some inconsistency. 

There have been several attempts to resolve the problem, by means of more sophisticated 
Ansatze for the full vertex |^-[10||, and by use of a so-called nonlocal gauge-function |i~8H20|, 
but none of them are completely satisfactorily. Also other methods, such as the inversion 
method plf , e-expansion p2fl , and lattice calculations ||, have not given a final answer to 
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this question, although the lattice results favor the existence of a finite critical number of 
fermion flavors l5,23||. The difficulty is that it is numerically very difficult to observe the 



exponential decrease of the dynamical mass for increasing N found in |6|rjJ . 

In our paper, the existence of a finite critical number of flavors is confirmed by analyzing 
the coupled Dyson-Schwinger equations, for both the fermion and the photon propagator, 
using several different approximations for the full vertex. We show that, within a certain 
class of vertex Ansatze, there is a critical number of fermion flavors for dynamical mass 
generation in QED3. This number does not depend very strongly on the precise form of 
the vertex Ansatz. Also the general behavior of the fermion propagator does not depend on 
the exact details of the Ansatz. The essential point is to take into account the full vacuum 
polarization in a consistent way. 

In the next section, we describe the formalism in more detail, together with our trunca- 
tion scheme. Before addressing the question of chiral symmetry breaking, we first analyze 
the full fermion propagator in the massless phase, see Sec. [TTI|. This already yields some non- 
trivial results, namely that the structure of the massless full fermion propagator is almost 
independent of our Ansatz. In Sec. |TV], we discuss dynamical mass generation, and show 
that there is a critical number of fermion flavors for chiral symmetry breaking, A c ~ 3.3, 
almost independent of our choice for the vertex. Finally we give some concluding remarks 
in Sec. |V|. 



II. FORMALISM 

A. Three-dimensional QED with N massless fermion flavors 

We consider QED3 with N massless fermion flavors, and choose to work in Euclidean 
space, ignoring the issues discussed in The Lagrangian in a general covariant gauge is 
given by 

N 

c = y: w 9 + e m + \f% + ■ (!) 

We use four-component spinors for the fermions, and accordingly a four- dimensional repre- 
sentation for the 7-matrices. With such a representation we can define chirality just as in 
four- dimensional QED. This chiral symmetry can be broken dynamically by generation of a 
mass for the fermions. In this formulation, there can also be a parity breaking mass term, 
which conserves the chiral symmetry, but it is known that such a mass is not generated 
dynamically |16| . 

We have N fermion flavors, and consider both the large N limit, as well as the quenched 
limit, A | 0. In the quenched limit, the mass scale is defined by the dimensionful coupling 
constant e 2 , and there are no free parameters. Outside the quenched limit, N, the number 
of fermion flavors is the only free parameter. When using the large N expansion P[, we keep 
the product Ne 2 finite, and it is most convenient to define the mass scale by 

a = -5-- (2) 
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B. Dyson Schwinger equations for the propagators 

The Dyson-Schwinger (DS) equation for the fermion propagator is given by 

S-\p) = S \p) -e 2 J -f^^S(k)r„(p, k)D, u (p - k) . (3) 

Since parity is not broken dynamically, we can decompose the fermion propagator into 

S- 1 (p) = A(p)p'+B(p), (4) 
and rewrite the DS equation into 

A(p) = l + - J lTr[^^S(k)r v (p : k)D^(p-k)} : (5) 
r d 3 k 

B{p) = e 2 j — - 3 \Tr[ 1 ,S(k)r u (k,p)D, l/ (p- k)} . (6) 

The problem in analyzing this equation is the full vertex and full photon propagator. For 
the photon propagator we also have a DS equation, namely 

D~l{q) = D \ u {q) - ^ j ^^(k)T v (k,p- k)S(p - k) , (7) 

without introducing new unknown functions. It is more convenient to write this last expres- 
sion in terms of the vacuum polarization tensor Ii^ u {q) 

/d k 
—^S{k)T v {k,p-k)S{p-k). (8) 

Since the longitudinal part of the photon propagator is not affected by the interactions, 
because of gauge invariance, we can write the full photon propagator in a general covariant 
gauge 



with the vacuum polarization 11(g) defined by 

IW?) = (<V - ^ U(q) . (10) 

The vacuum polarization tensor has an ultraviolet divergence, which can be removed by 
a gauge-invariant regularization scheme. However, this divergence is only present in the 
longitudinal part, so by contracting U^q) with 

(11) 



4 



we can project out the finite vacuum polarization 11(g) [pjj . So, the coupled DS equations 
for the photon and fermion propagator form a set of three coupled equations for three scalar 
functions, and the only unknown function is the full vertex function. Note that both the DS 
equation for the fermion propagator, and the one for the photon propagator, are exact. 

In principle, we could write down a DS equation for the full vertex function as well, 
but this will not lead to a closed set of equations: the DS equation for the vertex involves 
a four-point function, and so on. The full set of DS equations forms an infinite hierarchy 
of coupled integral equations for the Green's functions. In order to solve the DS equation 
for a particular Green's function, we have to truncate or approximate this infinite set of 
equations. For calculating the propagators, we must find a reasonable approximation for 
the full vertex function r M (p, k). 



C. Truncation scheme 

Now, what is a reasonable approximation for the full vertex function? The most simple 
and in some sense "natural" approximation is to take the leading order perturbative vertex 

(12) 

This truncation is commonly used in studies of the fermion DS equation, and it is usually 
referred to as ladder or rainbow approximation, since it generates rainbow diagrams in the 
fermion DS equation, and ladder diagrams in the Bethe-Salpeter equation for the fermion- 
anti-fermion bound state amplitude. In principle there is a systematic way to improve this 
truncation, namely by taking the next-to-leading order vertex function, and so on. 

The obvious disadvantage of this approach is that it relies on perturbation theory, 
whereas the DS equations are nonperturbative equations. Until one performs the com- 
plete next-to-leading order calculation, one does not know how reliable the leading order 
calculation is. Even after obtaining such a next-to-leading order, there remains some doubt 
about the validity of this approach, due to the non- uniformity of the 1/N expansion. 

Furthermore, as a consequence on the one hand of using a nonperturbative method to 
calculate some Green's functions, and on the other hand of employing perturbation theory 
for other Green's functions, one violates the Ward-Takahashi (WT) identities relating these 
Green's functions, and loses gauge covariance. To be specific, the WT identity relating the 
full vertex to the fermion propagator 

q fl T fl = S- 1 (p)-S- 1 (k) } (13) 

is an exact identity, which holds order-by-order in perturbation theory, but if one uses per- 
turbation theory to approximate the full vertex, and the truncated DS equation to calculate 
the fermion propagator, one will (in general) violate this identity. 

Instead of using a perturbative approximation for the vertex function, one can also use 
other nonperturbative information to make an Ansatz for the full vertex. The full vertex 
can be decomposed into 12 different Lorentz structures, and four of them are uniquely 
determined by the WT identity f24j[ . So by imposing the WT identity, we can write the full 
vertex as 
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r„ = \{A(p) + A(k))i„ + + K) Mp J, _ k) 

+ + (14) 

where the dots represent the part of the vertex not constrained by the WT identity. Now 
one can simply neglect that unconstrained part, and take the above expression as Ansatz 
(usually called Ball-Chiu vertex) for the full vertex. However, a perturbative calculation 
shows that the unconstrained part is not zero |25| , and by just neglecting it you never see 
what you are throwing away. 

The WT identity is not the only requirement one can impose on the full vertex function. 
Other requirements are that it reduces to the bare vertex in the weak coupling limit, and 
multiplicative renormalizability also restricts the full vertex, even in a superrenormalizable 
theory like QED3. Furthermore, any Ansatz should have the correct symmetry properties. 
However, all these constraints do not uniquely determine the full vertex, they all leave parts 
of the vertex undetermined. So the results might depend heavily on the particular choice 
for the vertex. Another disadvantage of using a vertex Ansatz is that it is not possible to 
improve the approximation in a systematic way. 



D. Vertex Ansatze 

In this paper, we use several different Ansatze for the full vertex, including a bare vertex, 
and compare the results to see how much influence these different choices have on the 
propagators. For the sake of simplicity, we consider the form 

r>, k) = f(A(p), A(k), A(p - fc))7„ , (15) 

with the restriction that for A(p) = 1, the full vertex reduces to the bare vertex. This 
automatically ensures that in the weak coupling limit, the full vertex reduces to the bare 
one. Such an Ansatz has the same tensor structure as the bare vertex, and by restricting 
ourselves to this tensor structure, we simplify the numerical integrations considerably. It is 
also generally expected that this tensor structure plays the most dominant role in the full 
vertex. 

To be more specific, we will analyze the coupled DS equations for the propagators, using 
the following choices for f(A(p), A(k), A(p — k))\ 

1. the bare vertex: f(A(p), A(k), A(q)) = 1 

2. a simple Ansatz inspired by the Ball-Chiu vertex Eq. ([14]) f24|: / = \{A{p) + A(k)) 

3. / = A(p)A(k)/A(q) 

4. / = \{A(p) + A{k)f 

5. / = A(p)A(k). 
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Of these choices, only the first two are actually motivated on physical arguments. The last 
one is primarily motivated by the fact that the resulting truncated DS equations are easy to 
solve analytically (at least in the symmetric phase). The choices |3] and ^ are merely added 
in order to see how much the results depend on our particular choice; Ansatz |] is based on 
|5], but linear in the wavefunction renormalization, as it is supposed to be, whereas Ansatz |] 
is the quadratic form of Ansatz [|. The last two Ansatze are also inspired by the suggestion 
in H that the effective vertex correction in the fermion DS equation is quadratic in the 
wavefunction renormalization. 

Finally, we have to choose a gauge, and since the Landau gauge is the most convenient 
and commonly used gauge, we use it. Using the Ansatz for the vertex we have just described, 
the DS equations for the propagators reduce to 

a/ \ _ -i | 2e 2 f d 3 fc A(k)(p-q)(k-q)f(A(p),A(k),A(q)) 

{P) p 2 J (2vr) 3 A 2 (k)k 2 + B(k) q 2 (q 2 + U(q)) ' 1 ' 

p M o 2 f d 3 fc B{k) f(A(p),A(k),A(q)) ' 

{P) J (2ir)3 A 2 (k)k 2 + B(k) q 2 + Tl(q) ' 1 ' 

n(a)-Nr*[^(2e-4k a ^ ' Q ) m A(p) f (A(p) , A(k) , A(g)) 

LL{q) - ive j (27r)3 yzK 4« • q f J + ^ ^ 2( ^ 2 + ^ , U»j 

with q = p — k. 



III. SYMMETRIC PHASE 

First, we consider the massless fermion phase, so B(p) = 0, which is always a solu- 
tion of Eq. (fL7|). This reduces the problem to solving two coupled equations, one for the 
wavefunction renormalization A(p), and one for the vacuum polarization H(p). 

A. Analytical results 

With B(p) =0, the vacuum polarization reduces to 

nf . _ N 2 f d 3 fc 2k 2 - 4k ■ q - Q(k ■ q)/q 2 f(A(k + q), A(k), A(q)) 
{q> 6 J (2tt) 3 k 2 (k + q) 2 A{k)A{k + q) ' 1 ' 

We first consider Ansatze |3| and ^ because the vacuum polarization can then be calculated 
analytically. With Ansatz |3|, the vacuum polarization becomes 

n(g) = Ne 2 q/(8A(q)) = aq/A(q) , (20) 

whereas Ansatz |] leads to 

U(q) = Ne 2 q/8 = aq . (21) 

In the massless phase, the equation for the wavefunction renormalization becomes 

a/ \ _ 1 | 2e 2 r d 3 fc f(A(p),A(k),A(g)) {p-g){k-q) 
[P) p 2 J (2tt) 3 A(k)k 2 q 2 (q 2 + 11(g)) ' 1 ' 
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For the Ansatze ^| and [5], this can now be reduced to 



A(p) - 1 + MP)^J ^ A(q)qi + aq iP fcV ~ ' (23) 

^^^wlwr^rW 1 ' (24) 

where we use a = Ne 2 /8 rather than e 2 to define the energy scale. This last equation can 
be solved exactly 

A -i, \ n 2a /- 00 dq ( 2 2 g 4 - p A p + q \ . . 

Np l -K l Jo q z + aq \ 2pq \p — q\ J 

whereas Eq. ( p3|) gives a nonlinear integral equation for A(p) 

A-U \ 1 2a f°° d 9 ( 2 2 9 4 -P 4 i P + ?\ / ofiA 

A (p) = 1 - tt~<2 2 / I? \ 2 i — U p — o — ln i 1 ' ( 26 ) 

JMp z 7t z Jo A{q)q z + aq \ Ipq \p — q\J 

and thus gives an implicit relation for the wavefunction renormalization. 

However, both Ansatze lead to the same explicit expression, if we make one further ap- 
proximation, which is commonly used in this context. In the infrared region, q <C a, the 
vacuum polarization dominates: the denominators, A(q)q 2 + aq in Eq. (p3|), and q 2 + aq in 
Eq. fl2"ip, both behave like aq, whereas for large momenta, q ^> a, the wavefunction renor- 
malization is almost equal to one, and the contribution to these denominators coming from 
the vacuum polarization can be neglected. So in both cases the wavefunction renormalization 
behaves like 

w(/;s+/:S(« 2 - 2 -^H) 

1 f a 2 3p A + 6p 2 a 2 — a 4 ^ a + p 4 ^ \a 2 — p 2 
Ntc 2 \ 3p 2 Qp 3 a \a — p\ 3 p 2 




B. Numerical results 

We can also solve the integral equations for the wavefunction renormalization and the 
vacuum polarization numerically. Starting with bare propagators to calculate the vacuum 
polarization, we solve Eq. (|16|) with this vacuum polarization. Using this solution for the 
wavefunction renormalization, we again calculate the vacuum polarization and repeat this 
procedure until the solutions for both H(p) and A(p) converge to a stable solution. In Fig. [I], 
we have plotted the analytic solutions for the wavefunction renormalization as described 
in the previous section, together with our numerical result. We see that they agree very 
well with each other qualitatively, and that in the ultraviolet region the A(p) is almost 
equal to one, as expected, whereas in the infrared region, it deviates considerably from the 
perturbative value. In particular, it is important to notice that for p I 0, the wavefunction 
renormalization vanishes, as also can be seen from Eq. (|27|). 



S 



For the other vertex Ansatze, we cannot solve the integral equations analytically. Nu- 
merically, we find a very similar behavior for the wavefunction renormalization for all five 
different Ansatze, as is shown in Fig. 0(a). The vacuum polarization also seems to be quite 
insensitive to the vertex Ansatz, as can be seen from Fig. 0(b). So it turns out that the 
exact form of the Ansatz is not very relevant for the behavior of the fermion propagator, 
provided that one uses a consistent approximation scheme: consider the coupled equations 
for the fermion and photon propagators, and use the same approximation for the full vertex 
in both the fermion and the photon DS equation. However, the deviation of the vacuum 
polarization from the perturbative leading-order behavior 

U{q)=aq, (28) 

does depend strongly on the Ansatz, see Fig. ||(c). 

It can be explained that the different vertex Ansatze lead to the same behavior for the 
wavefunction renormalization. In the symmetric phase, Eqs. (|i6l)-(|T8l) reduce to 

a/ \ _ i , 2e 2 r d 3 k f(A(p),A(k),A(g)) (p ■ q)(k ■ q) 

[P) p 2 J (2tt) 3 A(k)k 2 q 2 (q 2 + 11(g)) ' 1 ' 

ttm - at.2 f d 3 k (2k 2 - 4fc ■ q - 6(fc ■ q)/q 2 )f(A(p), A(k), A{q)) 
lV{q) ~ 7Ve J (2^ A(k)A{p) k 2 p 2 • m 

In the ultraviolet region, the propagators, and also the full vertices, reduce to the bare 
ones, so the crucial region is the infrared region. In the (far) infrared region, the vacuum 
polarization dominates in the denominator of the kernel of Eq. (p9|), q 2 + 11(g), so we can 
approximate this denominator by 11(g). We can now observe that there is a cancellation 
in the infrared region between the different effects of the vertex Ansatz: the function / 
enters the integral equation both in the numerator, and, via the vacuum polarization, in the 
denominator. The remaining equation for the wavefunction renormalization resembles the 
Eqs. fl23|) and fl2"4|), which we have analyzed in the previous subsection. This cancellation 
explains why we get a very similar behavior for the wavefunction renormalization, using 
quite different Ansatze for the vertex. 

Note that for such a cancellation, the vacuum polarization has to depend quite strongly 
on the precise Ansatz, and in Fig. 0(c) we see that the deviation from the perturbative 
behavior is indeed governed by the powers of A(p) in our vertex Ansatz. Using a bare 
vertex, the vacuum polarization becomes quadratic in 1/A, with a vertex linear in A, the 
vacuum polarization becomes linear in 1/A and with an Ansatz which is quadratic in the 
wavefunction renormalization, the vacuum polarization does not contain any powers of A, 
and is (almost) equal to the perturbative one. 

Finally, we show in Fig. || the wavefunction renormalization for several different numbers 
of fermion flavors. This shows that there is only a quantitative dependence on the number 
of flavors, as long as we stay in the symmetric phase. The other Ansatze yield a similar 
result. 



C. Infrared Behavior 

As can be seen from Eq. (p7|), the wavefunction renormalization vanishes in the infrared 
region, in contrast to what one might expect based on ordinary perturbation theory or the 
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1/N expansion. Such a vanishing behavior could very well be related to an anomalous 
dimension for the wavefunction renormalization as suggested in It has been argued 

that the naive perturbative behavior 



^pert(p) -l + T7T7-^ln (p/a) 



3iVvr 2 

is the first term of the build-up of an anomalous dimension 



with 



A(p) 



V 



a 
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3A% 2 ' 

An expansion around the origin of the analytic solution Eq. 



A(p) = l/ 1 



3A% 2 



In (p/a] 



(31) 
(32) 
(33) 

gives as leading contribution 

(34) 



So although we do not find the anomalous behavior Eq. fl32[) explicitly, our result is in 
agreement with it to leading order in 1/N. Moreover, our solution vanishes at p I 0, just 
like the suggested solution Eq. ([32]), and does not have the unphysical behavior of the 
perturbative result Eq. (|3~1~D, which diverges at small momenta. 

The results with the Ansatz confirm this anomalous behavior of the wavefunction 
renormalization very well. In Fig. |3|, we have also plotted the behavior Eq. ([32]), together 
with our numerical results using this Ansatz. This shows that our numerical results are 
in good agreement with the expectation of an anomalous dimension for the wavefunction 
renormalization, at least in the infrared region. In the ultraviolet region, Eq. (^) is not 
expected to be valid. 



IV. DYNAMICAL MASS GENERATION 

One of the interesting features of QED3 is that a fermion mass can be generated dynam- 
ically, breaking chiral symmetry. Starting with massless bare fermions, they can acquire a 
dynamical mass through nonperturbative effects. The order parameter for this symmetry 
breaking is the chiral condensate (ipip), but it is more convenient to consider the infrared 
value of B(p) as the order parameter. 

Writing the full fermion propagator as 

S-\p)=A(p)p'+B(p), (35) 

a nonzero solution for B(p) implies a nonzero condensate, and signals dynamical mass gen- 
eration. The infrared value of the dynamical mass function defined by m(p) = B(p)/A(p), 
m(0) = B(0)/A(0) can also be used as order parameter. Note we are not calculating the 
physical mass, defined at the pole of the propagator; this physical mass is expected to be 
in the Minkowski region, at least for observable particles, whereas we are performing our 
calculations completely in Euclidean region. The question of the existence of such a physical 
mass in the time-like region is not addressed here [13[ . 
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A. Existence of a critical number of flavors 



Based on bifurcation theory, we can show that, for determining the critical number (if 
there is a critical number) of fermion flavors, it is sufficient to keep the terms which are 
linear in the generated mass, assuming that there is a continuous phase transition. This 
means that we can use the symmetric solutions for the vacuum polarization and wavefunc- 
tion renormalization in the equation for B. To avoid infrared problems, we replace the 
denominator A 2 (k)k 2 + B 2 (k) by A 2 (k)k 2 + B 2 (0), and thus the equation for B reduces to[| 

p M „ 2 f d 3 k B{k) f(A(p),A(k),A(q)) 

{P) J (2ir) 3 A 2 (k)k 2 + B 2 (0) q 2 + Tl(q) ' 1 J 

or in terms of the mass function 

m(r>) = 2e 2 f d * k ^ f (Mp) , A(k) , A(q)) 

KP > J {2Tifk 2 + m 2 {Q)A(p)A(k)(q 2 + Ii(q))' 1 ' 

with A and II as found in the previous section. 

Consider first the Ansatze |3| and [5]. They lead to 

, . „ , f d 3 k m(k) 1 ,„ . 

m(p) = 2e 2 - — - — — ^. . . . „ , 38 

Ky ' J (2-KYk 2 + m 2 (Q)A(q)q 2 + aq 1 V ; 

and 

. . „ 9 f d 3 k m(k) 1 
m (p) = 2e 2 1 , 39 

J (2n) 4 k + m z {0) q z + aq 

respectively. We note that these equations are qualitatively similar: in the infrared region, 
the part coming from the vacuum polarization, aq, will dominate over the other part, A(q)q 2 
or q 2 respectively, whereas in the ultraviolet region A{q) ~ 1. Thus, both Ansatze lead to 
almost the same equation for the mass function. Note that Eq. ( P5|) is exactly the same 
equation for the mass function as one would get by using a bare vertex and neglecting the 
wavefunction renormalization completely [0]. 

Using the fact that the essential region for dynamical mass generation is the infrared 
region p C q, we can use the scale a as an ultraviolet cutoff. This reduces both Eq. (|38|) 
and Eq. (^) to exactly the same equation, namely 

8 r a dk k 2 m{k) 
m ^ V) = Wri 2 Jo max(p, k) k 2 + m 2 (0) ' ^ ' 

It is well-known that this equation leads to a critical number of fermion flavors for dynamical 
mass generation. This is seen directly by reducing it to a second-order differential equation 



1 Alternatively we could neglect the mass function in the denominator and introduce an infrared 
cutoff, which is generally identified with the infrared value of the mass function. 
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(41) 



with boundary conditions. Below the critical number of fermion flavors, N c = 32/n 2 ~ 3.24, 
the solution is given by a hypergeometric function 



P 2 



m(p) = m(0) 2 Fi [a+, a_, §; -§-^ ] , (42) 



m 2 (0) 

with 



a ± = -± iV /^c/iV-l. (43) 

As iV approaches the critical number of fermion flavors, the infrared value of the mass 
function m(0) decreases rapidly according to 



m(0) = a exp 



-2tt -, 
+ 31n2 + ^7r 



iV c /iV 



2 



(44) 



Above the critical number of flavors, the trivial solution m(p) = is the only solution. 

The equation for the dynamical mass function using the other Ansatze cannot so easily 
be obtained. However, using a simple counting argument, we can already expect that there 
is a critical number of fermion flavors, independent of the precise form of our Ansatz. In 
the infrared region, which is essential for dynamical symmetry breaking, we can neglect 
q 2 with respect to the vacuum polarization in the denominator q 2 + 11(g). Counting the 
powers of A, we see that this vacuum polarization is roughly proportional to / /A 2 , so the 
dependence of the integration kernel in Eq. ([37]) on the wavefunction renormalization and 
the vertex Ansatz might cancel out. This is similar to what we have already seen in the 
previous section, namely that the behavior of the wavefunction renormalization is almost 
independent of the vertex Ansatz we use. 

Therefore, although the behavior of the vacuum polarization might depend on the vertex 
Ansatz, we expect that in fact the fermion propagator is not very sensitive to the precise 
Ansatz, due to a cancellation in the (far) infrared region between the implicit dependence 
of n(g) on the function f(A(p), A(k), A(q)), and the explicit appearance of this function in 
the integration kernel of Eq. fl3"T|). On these grounds one might expect a critical number of 
fermion flavors, independent of our choice for the vertex. 

Even if we relax our requirement that the Ansatz goes to the bare vertex in the ultraviolet 
region, we find a similar result: consider for example the (unphysical) vertex 

2A(p)A(k)^. (45) 

This would lead to the equation for the mass function like 

2 f m(k) 2 

m{p) = 2e 7(2^ P + mW + 2cV ( 6) 

which has an additional factor of two. However, in the infrared region, this factor of two 
cancels against the additional factor of two which comes from the vacuum polarization, 
leading to the same critical number as with Ansatz |5|. 
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Note that it is crucial for such a cancellation to occur that we consider the vacuum 
polarization with full propagators; using bare propagators in the loop for the vacuum po- 
larization, is an inconsistent approximation. Formally, using such bare fermions in the loop 
for the vacuum polarization is in agreement with the 1/N expansion, but as we have seen 
in the previous section, the wavefunction renormalization is not of the order 1 + 0(1/ N) 
in the infrared region. Therefore we should use the full propagators when calculating the 
vacuum polarization. Even when using bifurcation theory to calculate the critical number 
of fermion flavors, one should use the full propagators in the vacuum polarization, or to be 
more precise, the massless full propagators. 



B. Numerical results 

We can solve the coupled equations numerically for the mass function m (or for B), 
the wavefunction renormalization A, and the vacuum polarization II in the broken phase. 
Starting with a trial function for B and the leading order contribution for the wavefunction, 
A(p) = 1, we can evaluate the vacuum polarization and solve the coupled equations for A 
and B. Next, we calculate the vacuum polarization, using these numerical solutions, and 
iterate this procedure until all three functions converge to a stable solution. 

For the quenched approximation, N = 0, all five Ansatze lead to the same solution for 
B, since the vacuum polarization is zero in this case, and the wavefunction renormalization 
equal to one. This mass function is almost constant in the infrared region, and decreases 
rapidly, as 1/p 2 , at large momenta, for p > e 2 . It agrees very well with earlier numerical 
studies of this approximation [|]. 

The typical behavior of the functions A, B, and II is shown in Fig. [|for several values of 
N. We note that the wavefunction renormalization decreases for small momenta, but does 
not vanish for p J, 0: its value at the origin is nonzero, in contrast to the behavior in the 
symmetric phase. Also the behavior of the vacuum polarization is different: in the broken 
phase, the vacuum polarization behaves like p 2 in the (far) infrared region, so U(p)/p 2 is 
finite at p I 0, whereas in the symmetric phase, U(p)/p 2 blows up at the origin. The infrared 
behavior of the vacuum polarization is governed by the generated mass function: the smaller 
this mass function, the larger U(p)/p 2 . Finally the mass function, or rather B, is almost 
constant at small momenta, and decreases rapidly for large momenta. However, in the 
unquenched approximation, and especially close to the critical number of fermion flavors, 
we see two relevant mass scales pD| : B(p) starts to decrease at the scale of the generated 



mass, but only beyond the energy scale a, it decreases like 1/p 2 . This phenomenon might 
be relevant for hierarchy problems in unified theories. 

To determine the critical number of fermion flavors, we consider N as a continuous 
parameter in the DS equations. In Fig. |5|, we show the infrared values of the wavefunction 
renormalization, of the scalar function B, and of the vacuum polarization, or rather of 
U(p)/p 2 , as functions of N. For increasing N, we note that B(0) decreases, and that the 
infrared value of U(p)/p 2 increases rapidly. Both functions indicate a critical value of the 
number of fermion flavors for dynamical mass generation, at which B vanishes, and where 
the limit p J, of U(p)/p 2 diverges. Also the wavefunction renormalization at the origin 
decreases quite rapidly close to this critical number of flavors. 
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The dependence of m(0) as function of iV at fixed a is shown in Fig. |6|. For comparison, 
we also included the analytical result Eq. (fHj). We see very clearly that the generated 
infrared mass also decreases very rapidly for increasing N, and the figure indicates a critical 
number of fermion flavors of about 3.3, more or less independent of the vertex Ansatz. The 
Ansatze |3], f|, and [5] give exactly the same critical behavior, and also Ansatz ^| leads to 
almost the same results. Only a bare vertex, Ansatz [I], indicates a slightly higher value for 
the critical number of flavors, but the general behavior is the same for all five Ansatze. 



V. CONCLUSION 

We have solved the coupled DS equations for the fermion propagator and the vacuum 
polarization, both in the chirally symmetric and in the broken phase, using a certain class 
of vertex functions. With this type of vertex, the behavior of the fermion propagator is 
almost independent of the exact form of the full vertex. We find a critical number for chiral 
symmetry breaking, N c = 3.3, below which there is dynamical mass generation; above this 
critical number, only the chirally symmetric solution exists. 

In the chirally symmetric phase, the wavefunction renormalization is approximately equal 
to one in the ultraviolet region, but in the infrared region, A(p) vanishes, indicating an 
anomalous dimension. This non-uniform behavior of the wavefunction renormalization lies 
at the origin of a controversy about the existence of a critical number of fermion flavors. 
Also in the broken phase, A(p) is not of the order one, as one might expect naively, but 
considerably smaller: as the number of fermion flavors approaches the critical number, A(0) 
tends to go to zero. These results will also be relevant for studies of this model at finite 



temperature [p7 |. 



Our main result, a finite critical number of fermion flavors for dynamical mass generation, 
confirms the assertions obtained earlier by Appelquist et al. @. Their analysis was based 
upon a bare vertex and neglect of the wavefunction renormalization. There have been 
several claims that including the effects of the wavefunction renormalization leads to a 
different result, namely dynamical mass generation for all numbers of fermion flavors ||[7|]- 
These and other studies, such as @,|], show a crucial dependence on the behavior of the 
full vertex; in general a bare vertex or a vertex linear in the wavefunction renormalization 
would lead to chiral symmetry breaking for all N, whereas a vertex quadratic in A leads 
to a critical number of fermion flavors. However, all these studies include the wavefunction 
renormalization and an Ansatz for the full vertex in the fermion DS equation, but not in 
the vacuum polarization. 

We would like to stress that in order to solve the DS equation for the fermion propagator 
self-consistently, we have to treat the fermion propagator nonperturbatively both in the 
fermion DS equation and in the equation for the vacuum polarization, at least if we are 
using the unquenched approximation. Also, one must use the same approximation for the 
full vertex in both the fermion and photon DS equation. It turns out that the effects of 
the wavefunction renormalization through the vacuum polarization change drastically the 
naive results obtained in In contrast to these earlier results, our results are almost 

independent of the vertex Ansatz. 

For studying the chiral phase transition, one can use bifurcation theory and therefore 
neglect the dynamical mass function in the wavefunction renormalization. However, one 
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should keep in mind that bifurcation theory does not imply that one can use bare propagators 
in the vacuum polarization: one should use the chirally-symmetric full propagators. Of 
course, if the wavefunction renormalization is equal to one (as one often assumes in this kind 
of calculations), one can use the bare propagators to calculate the vacuum polarization. 

We have made only one approximation in our calculation: replacing the full vertex by 
our Ansatz, Eq. (|T5|) . Of course, one could question this approximation, but our result is 
almost independent of the precise form of the function /. Furthermore, it has been shown 
that a next-to- leading order calculation in the context of the 1/N expansion also leads to a 
finite critical number of fermion flavors || . Our result is also in good agreement with lattice 
calculations ||. 
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FIGURES 



FIG. 1. The the numerical and analytical solutions for the wavefunction renormalization in the 
massless phase for N = 6, using the Ansatze || and |5[ 

FIG. 2. The numerical solutions in the massless phase for N = 6, using the Ansatze [l||5]: (a) the 
wavefunction renormalization, (b) the vacuum polarization, and (c) the deviation of the vacuum 
polarization from the perturbative behavior. 

FIG. 3. The wavefunction renormalization in the massless phase for N = 2, 4, 6, and 10, using 
Ansatz ||, together with the behavior based on the existence of an anomalous dimension. 

FIG. 4. The functions A(p), B(p), and Ii(p)/p 2 in the chirally broken phase for N = 0, 2, and 
3, using Ansatz |2|. 

FIG. 5. The infrared values of the functions ^4(0), B(0), and lim p ^o n(p)/p 2 as functions of N 
at fixed e 2 

FIG. 6. The infrared values of the mass function m(0) as function of iV at fixed a. 
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